arXiv:1504.07928v4 [hep-ph] 24 Sep 2015 


CERN-PH-TH-2015-091 

TTP15-018 

Lepton-flavour violating B decays in generic Z' models 

Andreas Crivellin,^ Lars Hofer,^ Joaquim Matias,^ Ulrich Nierste,^ Stefan Pokorski,"* and Janusz Rosiek^ 

^ CERN Theory Division, CH-1211 Geneva 23, Switzerland 
^ Universitat Autonoma de Barcelona, 08193 Bellaterra, Barcelona, Spain 
^Institut fiir Theoretische Teilchenphysik, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany 
^Institute of Theoretical Physics, Department of Physics, University of Warsaw 

LHCb has reported deviations from the Standard Model in 6 —> sfi^ transitions for which a 
new neutral gauge boson is a prime candidate for an explanation. As this gauge boson has to couple 
in a flavour non-universal way to muons and electrons in order to explain Rk, it is interesting to 
examine the possibility that also lepton flavour is violated, especially in the light of the CMS excess 
in /i —>■ . In this article, we investigate the perspectives to discover the lepton-flavour violating 

modes B —^ K^*^T^pA, Bs —^ jjA and B —^ K^*^, Bs —>■ For this purpose we consider 

a simplihed model in which new-physics effects originate from an additional neutral gauge boson 
{Z') with generic couplings to quarks and leptons. The constraints from r — >■ 3/r, t — >■ pLVV, p. — >■ ey, 

Qfj, — 2, semi-leptonic b —>■ sp^p~ decays, B —>■ K^*'‘uv and Bs-Bg mixing are examined. From 
these decays, we determine upper bounds on the decay rates of lepton flavour violating B decays. 

Br(B —>• Kvv) limits the branching ratios of LFV B decays to be smaller than 8 x 10“®(2 x 10“®) 
for vectorial (left-handed) lepton couplings. However, much stronger bounds can be obtained by a 
combined analysis of Bg-Bs, t —> 3p, t —> pvt? and other rare decays. The bounds depend on the 
amount of fine-tuning among the contributions to Bg-Bg mixing. Allowing for a fine-tuning at the 
percent level we find upper bounds of the order of 10“® for branching ratios into rp final states, 
while Bg —>■ p'^e^ is strongly suppressed and only B —> K^*^p^e^ can be experimentally accessible 
(with a branching ratio of order 10“^). 


I. INTRODUCTION 


While most flavour observables agree very well with 
their Standard-Model (SM) predictions, there are some 
exceptions in semi-leptonic B decays (see for example [T] 
for a recent review). LHCb [2] recently found indications 
for the violation of lepton-flavour universality in the ratio 


Br[B ^ Kp+p-] 
Bt[B Ke+e-] 


± 0.036, 


( 1 ) 


which deviates from the theoretically clean SM prediction 
= 1.0003 ± 0.0001 [3] by 2.6 cr. In addition, LHCb 
has reported deviations from the SM predictions mz] in 
the decay B K*p'^ p~ (mainly in an angular observ¬ 
able called P5 0) with a significance of about 3 a PEU]. 
Furthermore, also the measurement of Br[Hs —>■ cfp'^ p~] 
disagrees with the SM prediction [nUm by about 3cr p. 

Interestingly, these discrepancies can be explained in 
a model-independent approach by a rather large new- 
physics (NP) contribution Cg ^ to the Wilson coefficient 
of the operator (the component of the usual SM op¬ 
erator Og that couples to muons, see eq. ([|)) [MIS]- It 
is encouraging that the value for Cg ^ required to explain 
Rk (with Cf® = 0) is of the same order as the one needed 
for B —)• K*p'^p~ and Bg —t cfp'^p~ P ED]. Taking 
into account the 3fb“^ data for B —)■ K*p'^p~ recently 
released by the LHCb collaboration cni, the global sig¬ 
nificance is found to be 4.3 a for NP contributing to 
only, and 3.13cr in a scenario with Cg^ = —C^g [ID] . 

Many models proposed to explain the b —>■ sp~^p~ data 
contain a heavy neutral gauge boson (Z') which gener¬ 
ates a tree-level contribution to Cg^ [T3||2TH25]. If the 
Z' couples differently to muons and electrons, Rk can 


be explained simultaneously [25H29] . Since in this case 
lepton-flavour universality would be violated, it has been 
proposed to search for lepton-flavour violating (LFV) B 
decay modes as well EOj. This is also motivated by the 
CMS excess in Br[/i —>• pr] [3T| which can be explained 
simultaneously together with Rk, Br[Hs —>■ (j)p'^p~] and 
Br[H —>• K*p'^p~] within a single model [DDIET] . 

While the specific model of Refs. EillZ] predicts only 
small effects in LFV B decays, the situation could be dif¬ 
ferent in a generic model. In this article we examine the 
LFV decays B —>■ K^*'>T^p^, Bg -T p^ (and the cor¬ 
responding p^e^ channels) studying a simplified model 
in which the NP effects originate from a heavy new gauge 
boson Z' of mass Mz' with generic couplings to quarks 
and leptons [DDj. We introduce the relevant Z' couplings 
to sb and charged lepton pairs f, f! = t, p,e via 

Cz' A + T^bSU^PLb + L^ R. (2) 

As the Z' is assumed to be much heavier than the scale 
of electroweak symmetry breaking, its couplings must re¬ 
spect SU{2)l gauge invariance. This implies that the 
couplings to neutrinos and to left-handed charged lep¬ 
tons are equal: To study bounds on 

the LFV B decay modes, we perform the following steps: 

1) Motivated by the model-independent fits to H —> 
K*p^ p~, Bg —)■ (pp'^p~ and Rk we consider two scenar¬ 
ios for the Z' couplings to leptons: scenario 1 assumes 
vectorial couplings, i.e. r|^, = F^, = corresponding 
to Cfg = C’lg = 0. Scenario 2 considers left-handed 
couplings, i.e. F^, =0, corresponding to Cg^ = —C^q . 

2) We use the experimental upper bound on i? — 
decays to set upper bounds on LFV B decays, indepen¬ 
dently of the values of ■ 
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FIG. 1: Feynman diagrams illustrating the steps 1-4 of our 
analysis (see text). The diagrams display the dominant Z' 
contribution to Bs — Bs mixing, B —>■ K /j,'^fj, , Bs —>■ 
fi~, T —> 3/i, r —>■ nvv and B —>■ . 

3) From Bg-Bg mixing we obtain upper limits on F^j, as 
a function of a fine-tuning measure (to be defined later). 

4) In the lepton sector the Z' couplings can be con¬ 
strained by T —>■ 3/x and r —>■ ^vv. 

5) Taking into account the constraints 3) and 4) we de¬ 
rive upper limits on the branching ratios of Bg —>■ 

B which are stronger than the ones ob¬ 

tained in 2), but depend on the amount of fine-tuning in 
Bg-Bg mixing. 

In Fig. [^we show the Feynman diagrams for the dom¬ 
inant Z' contribution corresponding to the steps 1-5 of 
our analysis. We apply a similar procedure to fi¬ 

nal states. In this case the best bounds on the lepton 
couplings are coming from /r —>■ ey and /i —> evD. 


II. PROCESSES AND OBSERVABLES 

In the subsections A-E we collect the formulae for the 
steps 1-5 of our analysis outlined in the introduction. 


A. Bg — Bg mixing 

Using the notation of Refs. [3H [33] for the operators 
describing Bg — Bg mixing, the first diagram in Fig. 
feeds the Wilson coefficients of 

Ol = [Sa"f^PLba] [spl^PLbls] , 

O 5 = [SaPLbp] [SjsPRba] , (3) 

as well as 0 [ obtained from Oi by interchanging Pr o 
Pr. The coefficients are 

cV = (r^f ^)'/(2M|,), Cs = -2Ff,F«/(M|,). (4) 


For QCD renormalization group effects we use the next- 
to-leading order equations calculated in Refs. [321133] . 


B. b —>■ sl'^ P transitions 


For b —>■ s£~^P transitions we need the operators 

< 10 ) = ^[srPLb][h,h^)£'], (5) 


and their primed counterparts found by Pr Pr- Z' 
contributions to other operators (such as the magnetic 
operator O 7 ) are negligible. The diagrams of Fig. give 


r^'W 

'-^9,10 


TT 1 

y2M|, aGpVtbVt* 


pL(i?) 
^ sb 


,( 6 ) 


which have to be multiplied by —AGFVthVtsl^f^ f^e 
effective Hamiltonian. 

As first noted in Ref. [TS] [M] a good fit to R —>■ 
data, leaving Br[Rs —>• unchanged, is 

obtained with Cg^ < 0 and ~ 0. Another 

interesting solution is given by Cg^ = —C^q [HlfTHj. 

In our analysis we use the global fit of Ref. jb] [18], 
resulting for the two scenarios under consideration in 


- 0.53(-0.81) > > (-1.32) - 1.54, (7) 

-0.18 (-0.35) > = >(-0.71) -0.91,(8) 

at the (1 cr) 2 a level, respectively. The quoted ranges are 
in good agreement with preliminary results of Ref. m- 
Note thatBr[i?s —)■ ] is suppressed in scenario 2 

compared to the SM. This effect is taken into account 
via the global fit used in our analysis. 


C. B ^ 


Following |i 
nian as 

TTW' _ 

PeS — ~ 


we write the relevant effective Hamilto- 


4Gf 

71" 


Vtb, v^g 


-I 

47r 


/-iizv __ _ _ _ _ 

“ ^/2M}, aGpVtbVt^g 


?r'+c7o7) 

(9) 


(10) 

Y^L{R) -pij 

^ sb ^ vv' ’ 

(11) 


In the approximation F^j, = 0, the branching ratio (nor¬ 
malized to the SM prediction) reads 


R 


VIZ 

km 



( 12 ) 


with « — 1.47/s^ « —6.4. The complete expres¬ 
sions can be found in Ref. |35j . The current experimental 
limits are R'^ < 4.3 [33] and < 4.4 [37| . 

Due to SU{2) invariance, we have = {Cg^ — Cl{j)/2, 
so that Cfi = Cg^ j 2 in scenario 1 and = Cg^ in 
scenario 2 . 
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D. T —>■ ^vv, /X —>■ evu and t —>■ 3/x 

The Z' boson contributes to r —>■ in two ways: it 

generates loop corrections to the W exchange diagram (as 
in the lepton-flavour conserving case [25) 1 and it mediates 


T —>■ /xi/P at tree-level via LFV couplings. The latter 
contribution decouples as 1 /m^, from the branching ratio 
Br [r —>■ /xz/jx] for final-states where it interferes with 
the SM tree-level amplitude, and as 1/ m^, for other final- 
state flavours We find 


Br \t 


3r^„r:(:'^ logm^/m|, 




[t /xzxz/] =Br [t 1 H -^- /- t — T" 

^ ^ ^ JSM 1 4^2 1536 v^ 7 r 3 r^TO|, 


Re 


■pij 


+o 


■ (13) 


The HFAG value [38] for the branching ratio reads 

BR(t ^ /xzx^i>^)exp = (17.39 ± 0.04)% . (14) 

This should be compared to 

BR(r ^ piVrV^)sM = (17.29 ± 0.03)% , (15) 

obtained from the SM prediction in Ref. [39] and a com¬ 
bination of the T lifetime measurements in Refs. [101145) . 
The difference is given by 

= Br(r - Br(r ^ P-Vri'f,)exp 

= (-1.0 ± 1.1) X 10"% (16) 

at the 2 a level, adding the error originating from the SM 
theory predictions linear to the experimental one. In the 
analogous case of F^e we demand 

< 4 X 10-% (17) 

This choice restricts corrections to the Fermi-constant, 
defined through the decay /i —)■ euP, to the sub per-mille 
level and thereby avoids conflicts with electroweak preci¬ 
sion data. 

The Z' boson further mediates the LFV three body 
decay r —>■ 3/x at tree-level, with the branching ratio given 


by (cf. e.g. [46ll47] l 


Br [r —)• 3/i] = 


m" 


|-pi? pi? I 

I ^ fifi \ 


15367r3r^M|, L 

-b |f^ f^ I' 

' I f^T flfl I 




Ipi? pL I 

I ^ flT^ (1(1 \ 


(18) 


Combining Belle [48] and BaBar [49] data gives 
Br [r —)• 3/i] < 1.2 x lO”® at 90% C.L. [35]. The corre¬ 
sponding decay /i —>■ 3e does not affect our phenomenol¬ 
ogy, because it involves Fee which we set to zero to com¬ 
ply with Rk- 


E. Lepton-flavour violating B decays 


Here we give formulas for the branching ratios of LFV 
B decays, taking; into account the contributions from the 
operators and relevant for our model. For 

Bs — )■ (with I ^ H!) we use the results of Ref. [SO] 

neglecting the mass of the lighter lepton. The branching 
ratios for B , B —>■ are computed 

using form factors from Ref. m (see also Refs. mm)- 
The results read 


with 


Bt[B,^ £+£'-] = 


2^2 rr *,2 f, Ma.x[mj,m‘j: 

-a Gp iVtbVtsl 1 ~ 


647r3 


Ml 


^9 ^9 


^10 ^10 


Br[H ^ = 10-® {apMu, 


O 9 -h Og 




riti' I 

Olo -h Olo 


2 2 
I ^ I j /^li' 

"T Og — O 9 -h CLK(*)a' ^10 ~ ^10 



a-KU' 

bKU' 

CkU' 

d-KW 

a-K-U' 

bK*U' 

ck-w 

dK-W 

r/i 

9.6 ± 1.0 

10.0 ± 1.3 

0 

0 

3.0 ±0.8 

2.7 ±0.7 

16.4 ±2.1 

15.4 ± 1.9 

/xe 

15.4 ±3.1 

15.7T3.1 

0 

0 

5.6 ± 1.9 

5.6 ± 1.9 

29.1 ±4.9 

29.1 ±4.9 


(19) 


( 20 ) 
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Note that the results [64] in Eqs. (191 and (20) are for 
l ~final states and not for the sums l'^ + 

t^l'~ constrained experimentally [38] : 


Br [B+ 

Br [B+ 

Br [B 
Br [B, 


< 4.8 X 10"^ 

< 9.1 X 10"®, 

< 1.4 X 10"®, 

< 1.2 X 10"®. 


( 21 ) 


III. PHENOMENOLOGICAL ANALYSIS 


First of all, one can already derive an upper limit on 
LFV B decays from B ^ Kuu alone, simply by employ¬ 
ing gauge invariance [65] . As one can see from Fq. (12) 


the contribution for LFV couplings can only be positive. 
Therefore we can give a strict upper limit on \Cq'^ \ assum¬ 
ing that all other contributions vanish [66]. We obtain 
1^9^I < 46 for our scenario 1 and \Cq^\ = \Ciq\ < 23 for 
scenario 2. This results in upper limits on the branching 
ratios of & —?► st^ decays: 


Br[B K*Tfi] 


< 


Br[i3s —>■ Tfj] « 2Br[i? — Air/i] 

8 X 10"® in scenario 1, 

2 X 10"® in scenario 2. 


( 22 ) 


However, as we will show now, even stronger constraints 
can be obtained by employing the combined constraints 
from the other observables. Let us first examine the nu¬ 
merical impact of the leptonic constraints. As seen from 
Fig. § for our scenario 1 (vectorial couplings), r —> /rz/i/ 
rules out an explanation of via a non-vanishing 
(contrary to claims in Ref. [53] where r —> pLvv was 
not considered). The constraints from Z —>■ and 

Z —>■ T^/i^ as well as from neutrino-trident production 
(NTP) (see Ref. [51]) are irrelevant in the displayed T^^- 
r^T- region for the considered Z' masses (around 1 TeV 
and above). The situation is similar in scenario 2 (left- 
handed couplings). In this case the interference with the 
SM terms in is always destructive, albeit small. 

The most stringent constraints on the couplings 
stem from Bg — Bg mixing. Using the 95% CL results on 
AmB, by the UTfit collaboration [551157] [67] one obtains 


> 

D 

H 


5. 

I-H 



r;;^/(»iz'/TeV) 


FIG. 2: Allowed 2 a regions in the F))), — F)))- plane from r —>■ 
HW for = 0 (blue), VXr = —2 (yellow), F)),- = 2 (green), 
T —>■ 3/i (red) and (light grey) for rriz' = 1 TeV. The 
dependence of the bounds on the Z’ mass is only logarithmic. 
Although NP effects move to the right direction, it cannot 
be explained within our model and we do not impose it as a 
constraint later on in our analysis. 


the blue contour of Fig. ^ are weakened if T^^ and T^ 
have the same sign with7rff,| < |rfj or |^| > iTf-J, 
as a consequence of cancellations in eq. ( [M| . At the 
2(7 level, current b —>■ data requires a substan¬ 

tial non-zero contribution to Cg^, eliminating the op¬ 
tion |rfj > |rfj. Fig. I illustrates the combined con¬ 
straints from b —>■ sfi'^ fi~ data mm for different values 
of (scenario 1). In principle there is no upper limit 
on |r^j| as long as 6 —)■ s/i"/x" data permits small but 
non-vanishing contributions to the primed operators Cg 
and/or CJolSH]- Therefore we quantify the degree of can¬ 
cellation in Fq. ([24]) by the following fine-tuning measure: 




(ri)^ + (rfJ2 _ bBT^.T^, 


2aB. 


MI.ARb 


[(rf/)2 + (r«)2] - 1 , 


(26) 


- 0.10 < ARb, = Ams./Amgf - 1 < 0.23. (23) 

One can now derive limits on T^^ and T^ via the relation 
aRb. = ^ [(rf',)2 + (rf,)2 - feB.rf.rf,] . (24) 

The coefficients oBs, only exhibit a weak logarithmic 
dependence on Mz' (about 3% when varying Mz' from 
1 to 3 TeV) and we use the values at Mz' = 1 TeV: 

aBjMl, « 5700 TeV"2 , 5 b, « 8.8. (25) 


Restricting Ab, to an acceptable value limits the maxi¬ 
mal size |r^j|. As we are exclusively interested in scenar¬ 
ios with Cg iQ Cg^K), we neglect (F^)^ in Fq. (26) and 
express T^j, in terms of Ab, and ARb, as 

IrfbI/ Mz' = \/ Ai?B, (1 + Xb^)I (2aBj <cb,-\/1 -|- Ab, . 

Note that we take all couplings real to comply 
with CP data in Bg — Bg mixing. Using the maximal 
IARbJ allowed by Eq. (23), we find 


The bounds resulting from Eqs. (23) and ([24| (shown by 


CB, = max [VARB,/2aBj « 0.0045 TeV"^. (27) 
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FIG. 3: Allowed regions in the '^sb/Mz' - rfb/Mz' plane 
from Bs-Bs mixing (blue), and from the fit of 

Ref. [6] to B —> K*fi~, Bs —>■ and Rk, with = 

±1 (red), = ±0.5 (orange) and = ±0.3 (yellow). 
Note that the allowed regions with positive (negative) F^j, 
correspond to positive (negative) F^^. The bonnds are shown 
for rriz' = 1 TeV but their dependence on the Z’ mass is only 
logarithmic. 



FIG. 4: Maximal value of Br[B —> (red), Br[B —> 

(blue) and Br[B 3 —> (green) in scenario 1 as 

a function of Cg'^ for a fine-tuning of = 100 (solid lines) 
and Xbs =20 (dashed lines). The bounds are shown for 
mz' = 1 TeV but their dependence on the Z’ mass is only 
logarithmic. 


Combining the bound on and Eqs. (131,(18) we de¬ 
rive upper limits for the coefficient Cg^: 


2 bdTT^FT-CB 

1^9 I - 5 4^4 IT/ w* 


max{Br[T ^ 3p]^^p} x 


\crf < A 


gev^TT^r^cl 


Bs 


maxjA^^^^p} X (1 ± . 


(28) 


(29) 


For scenario 1 we obtain = 16 and = 4, while 

for scenario 2 we get = 3 and = 1. 

The bounds from r —> pi/p only depend on the fine- 
tuning measure Xb^-, while those from r —> 3p also de¬ 
pend on the value of Cg^ (and C^q in scenario 2) de¬ 
termined from the fit to 5 —>■ sfj,'^fi~ data. The latter 
bounds disappear in the limit Cg^ —>■ 0, as in this case 
the Z'pp couplings may vanish so that the t ^ 3 fj, decay 
does not receive contributions from Z' exchange. 

From the upper bounds on Cg^g, we can finally deter¬ 
mine the maximal branching ratios for the LFV B decays 
with rp final states. They are shown in Fig. for sce¬ 
nario 1 with Xb^ = 20 and Xb^ = 100 (in scenario 2 they 
are a factor of 1/2 smaller). The kink in the curves oc¬ 
curs at the point where the Cg^g-independent constraint 
from T —>■ fivD becomes stronger than the constraint from 
T —^ 3fj,. One should note that the bounds presented in 
Fig.i which are given for mz' = 1 TeV, have only a 


weak logarithmic dependence on the Z' mass. 

Comparing these results to the experimental upper 
limits in Fq. (21), we see that the current experimen¬ 


tal sensitivity is still two orders of magnitude weaker. 
However, LHCb will be able to achieve significant im¬ 
provements in these channels. 

In the case of pe final states, the stringent bound from 
Br[p —ey] renders LFV B decays unobservable in the 
Gg^ region favored by current b —>■ sp+p“ data. For 
Cg^ —>■ 0, Br[il —can become relevant with 
its maximal size being constrained to 0(10“^) from p —>■ 
evv. 


IV. CONCLUSIONS 

In this article we have investigated the possible size of 
the branching ratios of lepton-fiavour violating B decays 
Bg —> fjA, Bg —>■ p^e^, B —>■ and B —>■ 

kG) e^fjA in generic Z' models. Motivated by the model- 
independent fit to 6 —>■ s transitions, we have considered 
two scenarios, one with vectorial (scenario 1) and another 
one with purely left-handed couplings (scenario 2) of the 
Z' to leptons. 

From Br(B —>■ Ki^i/) one obtains limits on the branch¬ 
ing ratios of LFV B decays of 8(2) x 10“® for scenario 
1(2) simply by using gauge invariance. However, even 
stronger bounds can be obtained by combining the lep- 
tonic constraints with a limit on the amount of fine tuning 
in the Bg—Bg system. For a fine-tuning of Xb^ ^ 100, we 
have found that still sizeable branching ratios of 0(10“®) 
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are possible in both scenarios for r/r final states, while for 
/xe final states they can only reach 0(10“^) in a region 
of parameter space disfavoured by the current data on 
B —>■ Bs —t and Rk- 

Note added: During the publication process of this 
article new LHCb results on Bs —>■ were re¬ 

leased, increasing the discrepancy compared to the SM 
to 3.5a 
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